package com.mvp.lt.airlineview.utils;

import com.amap.api.services.core.LatLonPoint;

import java.util.HashMap;

/**
 *
 * @author Administrator
 * @date 2017/9/8
 */
public class LatLngUtils {
    //圆周率 GCJ_02_To_WGS_84
    static double PI = 3.14159265358979324;

    /**
     * @return 转换为真实GPS坐标后的经纬度
     * GCJ_02转WGS_84
     * @throws <异常类型> {@inheritDoc} 异常描述
     * @author 作者:
     * 方法描述:方法可以将高德地图SDK获取到的GPS经纬度转换为真实的经纬度，可以用于解决安卓系统使用高德SDK获取经纬度的转换问题。
     */
    public static HashMap<String, Double> delta(double lat, double lon) {
        double a = 6378245.0;//克拉索夫斯基椭球参数长半轴a
        double ee = 0.00669342162296594323;//克拉索夫斯基椭球参数第一偏心率平方
        double dLat = transformLat(lon - 105.0, lat - 35.0);
        double dLon = transformLon(lon - 105.0, lat - 35.0);
        double radLat = lat / 180.0 * PI;
        double magic = Math.sin(radLat);
        magic = 1 - ee * magic * magic;
        double sqrtMagic = Math.sqrt(magic);
        dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * PI);
        dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * PI);

        HashMap<String, Double> hm = new HashMap<String, Double>();
        hm.put("lat", lat - dLat);
        hm.put("lon", lon - dLon);

        return hm;
    }

    //转换经度
    private static double transformLon(double x, double y) {
        double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(x * PI) + 40.0 * Math.sin(x / 3.0 * PI)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(x / 12.0 * PI) + 300.0 * Math.sin(x / 30.0 * PI)) * 2.0 / 3.0;
        return ret;
    }

    //转换纬度
    private static double transformLat(double x, double y) {
        double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(y * PI) + 40.0 * Math.sin(y / 3.0 * PI)) * 2.0 / 3.0;
        ret += (160.0 * Math.sin(y / 12.0 * PI) + 320 * Math.sin(y * PI / 30.0)) * 2.0 / 3.0;
        return ret;
    }


    private final static double a=6378245.0;
    private final static double pi=3.14159265358979324;
    private final static double ee=0.00669342162296594626;

    //gcj-02  to  wgs-84
    public static LatLonPoint toWGS84Point(double latitude, double longitude){
        LatLonPoint dev=calDev(latitude, longitude);
        double retLat = latitude-dev.getLatitude();
        double retLon=longitude-dev.getLongitude();
        dev=calDev(retLat, retLon);
        retLat=latitude-dev.getLatitude();
        retLon=longitude-dev.getLongitude();

        return new LatLonPoint(retLat, retLon);

    }

    //wsg84 to  gcj02
    public static LatLonPoint toGCJ02Piont(double latitude, double longitude){
        LatLonPoint dev=calDev(latitude, longitude);
        double retLat = latitude-dev.getLatitude();
        double retLon=longitude-dev.getLongitude();
        return new LatLonPoint(retLat, retLon);

    }


    private static LatLonPoint calDev(double wgLat, double wgLon){
        if(isOutofChina(wgLat,wgLon)){
            return new LatLonPoint(0,0);
        }
        double dLat=calLat(wgLon-105.0,wgLat-35.0);
        double dLon=calLon(wgLon-105.0, wgLat-35.0);
        double radLat=wgLat/180.0*pi;
        double magic=Math.sin(radLat);
        magic=1-ee*magic*magic;
        double sqrtMagic=Math.sqrt(magic);
        dLat=(dLat*180.0)/((a*(1-ee))/(magic*sqrtMagic)*pi);
        dLon=(dLon*180.0)/(a/sqrtMagic*Math.cos(radLat)*pi);
        return new LatLonPoint(dLat,dLon);
    }

    private static double calLat(double x, double y) {
        double ret=-100.0+2.0*x+3.0*y+0.2*y*y+0.1*x*y+0.2*Math.sqrt(Math.abs(x));
        ret +=(20.0*Math.sin(6.0*x*pi)+20.0*Math.sin(2.0*x*pi))*2.0/3.0;
        ret +=(20.0*Math.sin(y*pi)+40.0*Math.sin(y/3.0*pi))*2.0/3.0;
        ret +=(160.0*Math.sin(y/12.0*pi)+320*Math.sin(y*pi/30.0))*2.0/3.0;

        return ret;
    }

    private static double calLon(double x,double y){
        double ret=300.0+x+2.0*y+0.1*x*x+0.1*x*y+0.1*Math.sqrt(Math.abs(x));
        ret +=(20.0*Math.sin(6.0*x*pi)+20.0*Math.sin(2.0*x*pi))*2.0/3.0;
        ret +=(20.0*Math.sin(x*pi)+40.0*Math.sin(x/3.0*pi))*2.0/3.0;
        ret +=(150.0*Math.sin(x/12.0*pi)+300.0*Math.sin(x/30.0*pi))*2.0/3.0;;

        return ret;

    }


    private static boolean isOutofChina(double lat, double lon) {
        if(lon<72.004 || lon>137.8347){
            return true;
        }
        if(lat<0.8293 || lat>55.8271){
            return true;
        }

        return false;
    }
}
